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Abstract: In this paper, PZT actuators are used to study the dynamic behavior of a rotating 
disk structure due to rotor-stator interaction excitation. The disk is studied with two 
different surrounding fluids — air and water. The study has been performed analytically and 
validated experimentally. For the theoretical analysis, the natural frequencies and the 
associated mode shapes of the rotating disk in air and water are obtained with the 
Kirchhoff-Love thin plate theory coupled with the interaction with the surrounding fluid. 
A model for the Rotor Stator Interaction that occurs in many rotating disk-like parts of 
turbomachinery such as compressors, hydraulic runners or alternators is presented. The 
dynamic behavior of the rotating disk due to this excitation is deduced. For the 
experimental analysis a test rig has been developed. It consists of a stainless steel disk 
(r = 198 mm and h = 8 mm) connected to a variable speed motor. Excitation and response 
are measured from the rotating system. For the rotating excitation four piezoelectric 
patches have been used. Calibrating the piezoelectric patches in amplitude and phase, 
different rotating excitation patterns are applied on the rotating disk in air and in water. 
Results show the feasibility of using PZT to control the response of the disk due to a 
rotor-stator interaction. 
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1. Introduction 

The transverse vibration of rotating disks has been a relevant research topic for decades because the 
disks have various applications in engineering. Rotating disks are a common element in rotating 
machines such as circular saws, wafer cutting machines, disk brakes, grinding wheels, data storage 
disks and turbine rotors. One of the cases of interest involving disk-like structures are the 
turbomachinery components with small gaps between rotating and static parts. In this case the 
perturbations originated by the static parts (guide vanes) superposed with the perturbations originated 
in the rotating parts (rotating blades) lead to a pressure pulsation known as rotor-stator interaction or 
RSI. RSI excites the rotating and the stationary part of the machine with a superposition of excitation 
patterns that depend on the number of harmonic, guide vanes and rotating blades and with a frequency 
content that depends also on the rotating speed of the rotor. RSI has been reported as the main cause of 
a critical failure when the rotating disk- like structure was rotating in air [1,2] and when it was rotating 
in water [3]. To avoid large vibration problems caused by the RSI, it is of paramount importance to 
determine the dynamic behavior of rotating disk-like parts, which involves the study of the natural 
frequencies and mode shapes of the structure and the study of the excitation characteristic. 

Early work on the vibration of rotating disk was carried out by Campbell [4]. In this study the term 
"critical speed" at which the stationary wave appears on the disk was coined. The effect of rotation on 
the natural frequencies of the disk was introduced by Lamb and Southwell [5]. Their study was 
focused on a disk which rotated about its center with constant angular velocity. In [6] Southwell 
studied the vibration of circular disks clamped at their center. He considered the effect of a shaft that 
clamps the disk at its center, on the natural frequencies of the rotating disk. This effect was not 
considered in the previous study. Tobias and Arnold [7] used the linear plate theory to provide further 
explanation of the "stationary wave" effect. They concluded that the introduction of an imperfection 
reduced the magnitude of disk vibration under resonance conditions. Experimental studies with 
rotating disks have been carried out recently. Medigholi [8] showed the appearance of two resonances 
in the rotating disk (viewed from the rotating system) for each natural frequency of the disk stationary 
wave. In this case the disk was excited from the stationary frame. Ahn [9] could determine the mode 
shape of the rotating disk using two stationary sensors. He measured the phase difference between the 
signals of two proximity probes for this purpose. In this case, an electromagnet placed on the 
stationary system was the disk exciter. Some studies concerning rotating disks with different boundary 
conditions have been carried out by means of numerical simulation [10-13]. In [10] Heo studied the 
effect of the angular misalignment in the natural frequencies of the disk. Bauer [11] studied the effect 
of the attachment of the disk on the stationary part. Imperfections in bladed disks were considered by 
Pust [12,13]. Most of the mentioned cases deal with very thin disks that rotate in air or vacuum, but 
none of them studied the effect of rotation when the disk rotates with surrounding heavy fluids such as 
water. Kubota [14] studied this effect but only from the analytical point of view, so no experimental 
results were given. 
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The RSI phenomenon was described in [15,16], as a superposition of wake and potential effects. 
The combination of blades in the rotating part and guide vanes determine the excitation shapes that 
are applied in the machine [14]. Some recent studies have been developed studying this 
phenomenon [17,18], nevertheless all of them are centered in the study of the flow characteristics and 
not on the dynamic behavior of the structure itself. To study this effect from the rotating frame it is 
advantageous to use an excitation and measurement system placed on the rotating system (rotating 
disk- like part of the machine), since stationary sensors and actuators could affect the flow 
characteristics around the rotating part of the machine. 

To excite the disk with a rotating excitation, light exciters (that do not affect the mass of the disk) 
have to be attached on the rotating frame. Because electromagnets or shakers are very heavy and may 
affect the mass of the structure, light and thin piezoelectric patches can be used in this case. PZTs are 
used in many cases as exciters [19-24]. Yang [19] studied the governing equation of an elastic plate 
due to the excitation of one PZT. He also modeled the dynamic behavior of the actuator. In [20], 
Gomis presented a control law for a piezoelectric actuator considering the hysteresis. Some studies 
have also been found with more than one acting piezoelectric patch. Cheng [21] placed several patches 
on a plate and studied the effect of the added mass of the actuators. In [22] Sekouri used piezoelectric 
patches to excite a thin circular plate. Finally, Wang [23] and Yan [24] studied the feasibility of 
suppressing the aerodynamic flutter of a rotating disk. In both studies the disk was rotating but the 
actuators were placed on the stationary frame. Although in some of the mentioned studies PZT 
actuators are used to excite a rotating disk, the exciters are placed in the stationary frame and no 
studies have been found with multiple PZTs actuators placed on the rotating structure acting as modal 
exciters. Furthermore, in the mentioned studies PZT actuators are used to excite very thin rotating 
disks (thicknesses of less than 1 mm), which is in the range of rotating disks such as CD drives, DVD 
drives and other data storage disks, and no studies have been found with PZTs actuators exciting thick 
and submerged disks in water. 

In the present paper, the feasibility of studying the dynamic behavior of a rotating disk structure due 
to a rotor-stator interaction, by means of PZTs is analyzed. The main objective is to determine the 
dynamic response of a disk, rotating inside a casing with air or full of water, when it is excited with 
different excitation patterns that simulate the RSI. The remaining sections of this paper are organized 
as follows: Section 2 presents an analytical model that uses the thin plate theory for an annular plate, 
coupled with the pressure that the surrounding fluid exerts on the rotating disk to determine the natural 
frequencies and mode shapes and a one-dimensional modeling of the RSI. Section 3 describes the 
experimental setup used and the procedure to study this topic with PZT actuators. Section 4 includes 
numerical results for the amplitude of resonance of the first several natural frequencies of the rotating 
disk in air and in water and comparison between analytical and experimental method. Concluding 
remarks are given in Section 5. 

2. Analytical Model7 

Prior to the experimental tests with PZTs actuators exciting a rotating disk, an analytical model to 
predict the response of the rotating disk (in air and water) due to different rotating excitation patterns 
(that simulate the RSI) under resonance conditions, is presented. 
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2.1. Structural Response of a Rotating Disk in Air and in Water 



To study the natural frequencies and mode shapes of the rotating disk in air and in water, the model 
developed by Kubota [14] is used; particularized for a totally confined disk (Figure 1). The governing 
equation of the transverse vibration of a rotating disk in vacuum in cylindrical coordinates is: 

2 

w = 0 (1) 



d 2 w 
PdK tt-v + D 
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p D and h D are the density and thickness of the disk respectively, w is the displacement of the disk in 

Eh? 

n is the bending stiffness of the disk, where E and v the Young and 



the axial direction. D — 

12(l-u 2 ) 

Poisson modulus of the material. 



Figure 1. Model for a confined disk. 




For certain ratios r in /r out the first several natural frequencies correspond to mode shapes with no 
diametrical circles [25]. Considering only these mode shapes, which are also the more apt to be excited 
by the RSI [14], Equation (1) can be simplified for an annular disk, assuming that the vibration of the 
disk is uniform in the radial direction. In this case the vibration of the disk is studied in an averaged 
radius r 0 = ^Jr in ■ r out . The surrounding fluid is considered taking into account the pressure that the 

fluid is applying on the disk surface at r 0 [14]: 

d 2 w D*d 4 w 

In this Equation, D* is a parameter that depends on structural properties of the disk. Furthermore, 
for the considered mode shapes, the disk vibration at r 0 can be expressed as: 

oo 

w = V A n e jne e jMnt (3) 

n=±2 

where Inl indicates the number of nodal diameters in the mode shape and the sign of n indicates the 
direction of the travelling wave that appears on the disk. a> n is the corresponding natural frequency. 
In [14], Equations (2) and (3) are supposed to be valid also for the modes |n| = 1. Nevertheless, the 
discrepancy between experimental and analytical values for these modes shown in that study, suggest 
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to include only the modes with |n| > 1, where the analytical model predicts with good accuracy the 
experimental values. To calculate Pr 0 , it is assumed that the fluid on the tank is inviscid and 
incompressible. Therefore, a potential flow on the tank is considered. The solution of the potential 
function for each of the two fields in Figure 1 is given in [14]. Applying the Equation of Energy in the 
non stationary form [26], the term Vr 0 in the disk frame can be expressed as: 

Pr 0 = P F r 0 £ ^ • e jn9 • e^ 1 K + nn up ) 2 coth (^) + {<o n + nn down ) 2 coth (*^f™) (4 ) 

n=+l V 0 0 / 

p F is the density of the fluid and Q. up & (} down is the rotating speed of the flow with respect to the 

disk (Figure 1). These two parameters depend on the viscosity of the fluid, on the geometry of the tank 
that contains the disk and on the speed of the disk. The dependence of n rot with fi up and n down is 
approximately linear [26]. Since O fluid up & ^n u id,down are very complex to obtain analytically 
(the calculus involves Navier-Stokes equations) a CFD calculation can be performed to obtain these 
two parameters. 

When the surrounding fluid is air, since the density of the flow is very small, the term p To can be 
neglected in comparison with the terms of Equation (2). Therefore, the natural frequencies for this case 
can be obtained as: 

2 _ n * D * 

In [25] an alternative formulation to calculate the natural frequencies of an annular disk is proposed. 
Combining Equation (5) and the method proposed in [25] the value of D* can be obtained. 

According to Equation (5), the natural frequencies of the rotating disk do not depend on the rotating 
speed, for slow velocities. In other studies [11], the dependence of the natural frequency on the rotating 
speed of the disk, viewed from the rotating system is studied. For slow velocities it is shown that only 
a very slight increase is observed. From Equation (5), it can be also seen, that a value of n-positive and 
n-negative will give the same positive solution of w^ air . If the positive solution of n-positive and the 
positive solution of n-negative are introduced in Equation (3), it can be seen that the resulting mode 
shape, is a standing wave on the disk with n nodal diameters (Figure 2). 



Figure 2. n = +2, n = +3, n = +4 for the disk rotating in air, viewed from the rotating frame. 




For the case that the surrounding fluid is water, Equations (3) and (4) are introduced in Equation (2), 
and the resulting Equation is obtained as: 
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From Equation (6), it can be seen that when the disk rotates, fl up &fi down are different than zero 
and this causes that the positive solution of (^ n ,water f° r n-positive is different than the positive 
solution of (t) nwater for n-negative. This difference increases when fl rot is increased (since 
S2 up &/2 down also increases [26]). Substituting the two solutions obtained for each n, in Equation (3), it 
is observed that two travelling waves travelling in opposite direction at different frequencies appear. 
The first one (n-positive and lower frequency) travels in the same direction as the rotating speed of the 
disk and the second one (n-negative and higher frequency) travels in the opposite direction [14]. 

A physical explanation for this effect is the influence of the added mass of the fluid on a forward 
wave and on a backward wave. According to [14], the free vibration of an annular disk is the 
superposition of a forward and a backward wave, for each diametrical mode n. For the annular disk 
with steady surrounding fluid, the added mass effect of this fluid on the forward and on the backward 
wave is the same and therefore both waves will have the same natural frequency and the corresponding 
mode shape at this frequency will be the superposition of both waves, which is a standing wave. With 
a relative rotation of the surrounding fluid with respect to the disk, the added mass effect will be 
different for the forward than for the backward wave, since the relative velocity of the fluid with 
respect to the wave will be different depending on the rotating direction of the wave. This causes, that 
the frequency of the backward wave will be different than the frequency of the forward wave. In this 
case, for each diametrical mode n a pair of natural frequencies, which correspond to the forward and to 
the backward wave, will appear on the disk. A similar effect is shown in [27], for a fluid-conveying 
pipe with periodic boundary conditions. 

Increasing 0. disk will increase (2 up & fi down and this will enhance the mentioned effect, which 

means to increase the difference between both natural frequencies. For higher values of O disk than 
considered in this paper, some terms may be included in Equation (1) [11] (due to centrifugal and 
Coriolis forces) and therefore the analytical solution would be modified. Furthermore, higher velocities 
of the disk leads to low pressure areas what could generate vapor bubbles (cavitation) [28,29], 
changing the added mass effect depending on the type and amount of cavitation [30]. Therefore, with 
the presence of cavitation, the solution of Equation (2) becomes much more complex. 

The influence of the radial gap in the axial vibration is not considered in Equation (6). Nevertheless 
in [31] is shown that for short radial distances from the disk to the tank, the natural frequencies are 
decreased. In Equation (6), this effect can be considered increasing slightly r 0 , matching experimental 
and analytical results for the non rotating disk. 
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2.2. Rotor-Stator Interaction 

An unidimensional model for the Rotor-Stator Interaction is presented in [17]. When the rotating 
blades of the rotor pass in front of the static vanes of the stator (Figure 3) the pressure field in the gap 
between blades and vanes can be described as the superposition of all the combinations l,k: 

Pik( 0 > 0 = A ik ■ cos(lZ 0 B s + 0j) • cos{kZ b B r + 0 k ) for I = 1,2, ... , oo k = 1,2, ... , oo (7) 

Figure 3. Rotor-stator interaction. 




Pu 



This pressure field can be viewed from the rotating frame or from the stationary frame. In this 
study, since the structural response is viewed from the rotating frame, the excitation will 
be also considered from the rotating frame. Transforming Equation (7) in the rotating coordinate 
(9 r = 6 S — fit) this pressure field can be expressed as: 

A lk 

Pi.k = -^-cos(ZZ o n rot t - {kZ b - lZ 0 )9 r + 0i- 0 k ) 

+ -^cos(ZZ o n rot t - (-ZZ 0 - kZ b )B r + 0j + 0 k ) (8) 

Aik Aik 
= — cos(ZZ o n rot t - Yl B r + 0j - 0 k ) + — cos(ZZ o n rot t - Y2 O r + 0i + 0 k ) 

From Equation (8) it can be seen, that the excitation shape depends on the number of guide vanes 
and rotating blades (fi = kZ b — IZ 0 ) and (y 2 = —IZ 0 — kZ b ). This number represents the number of 
maximums and minimums in the pressure pulsation in a circle (Figure 4) and the sign of y x and 
y 2 indicates the rotating direction of the excitation. If it is positive, it indicates that the excitation 
rotates in the same direction (faster) than the rotating disk-like part. If it is negative it rotates in the 
opposite direction (slower). 

Usually the term of the pressure pulsation that contains y 2 is not relevant for the RSI, since the 
several first mode shapes of a disk like structure have a small number of diametrical modes [17]. 

The excited frequency depends only on the number of guide vanes and on the rotating speed of the 
machine (lZ 0 I2 rot ). For higher number of harmonics (l,k) lower amplitudes A lk are expected. 
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Figure 4. y = +2 pressure pulsations for a rotating disk. 




Rotating turbomachinery components that suffer the RSI, are designed to avoid resonances during 
their steady state. Nevertheless, during the acceleration or deceleration of the rotor, since the rotating 
speed changes, a resonance can occur if the natural frequency coincides with the excited frequency and 
the excitation shape with the mode shape. 

2.3. Forced Response of the Disk Due to Rotor Stator Interaction 



The frequency response function (FRF) is defined as the relationship between displacement at point 
p when a force is applied at point q for an arbitrary frequency co. Considering the response in 
resonance of the mode n, i.e., co = co n : 



li 



h 
h p i 



If only the response of the point 1 is studied: 



X 



1,71 



li 



'-Iq 



h 



P<7. 



(4 

{F} n 



(9) 



Sin' 



(10) 



2.3.1. Air 



When the disk is rotating in air at lower speeds, the mode shape is a standing wave on the disk 
(Equation (5) substituted in Equation (3)). It is considered, that the disk is excited in q equidistant 
points with respect to the origin of angles (point 1) and that the RSI pattern is equal in magnitude for 
all the q points (Equation (8)). The vector \h\x "' h lq ] n , is obtained for each of the considered 
natural frequencies co n , using the information of the mode shape (in this case a standing wave). For the 
excitation {F} n , the pattern deduced in Equation (8) is introduced. Since the excitation changes its 
phase for each position (not all the points moving in phase or in counterphase), complex numbers have 
to be used to calculate the response [32]. For each of the considered natural frequencies: 
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/ / / (g - 1) ■ 2jt \\ . . / / (g - 1) ■ 27T \\V 

■ ( cos I — 5 — j J - J ■ sin I — 5 — J J J 

0 n is an arbitrary complex constant that depends on the mode shape considered and j is the complex 
unity. From Equation (11), it can be deduced that to excite the structural mode +n, the relationship 
between number of exciters q, excitation mode y x has to be the following: 

n = ±y 1 ± Aq (12) 

In Equation (12), n,y lt q are entire and positive numbers that are defined by structural mode, 
excitation shape and number of equidistant exciters respectively and A is an arbitrary entire number 
(including 0). From this equation it can be deduced, that the structural mode ±n is excited with one 
exciter whatever the excitation shape is, as for one exciter this is not defined. When y x = ±n the 
structural mode is excited with any number of exciters q. Otherwise, when y^±n, only for the number 
of exciters that accomplish Equation (12) the structural mode +n is excited. 

2.3.2. Water 

When the disk is rotating in water at lower speeds, two travelling wave appears for each n 
(Equation (6) substituted in Equation (3)). The same assumptions made for the excitation of the 
rotating disk in air are assumed now. The main difference in this case, is that the structural mode shape 
is a travelling wave and to express [/in ■•• h lq ] n also complex numbers have to be used, since 
generally all the points are not moving in phase or in counterphase to each other. For each of the 
studied mode shapes, when lZ 0 D. rot = ro n : 
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In this case, to excite the structural mode +n, the following equation has to be accomplished: 

n = y 1 + Xq (14) 

As for the case that the disk rotates in air, for one exciter the structural mode +n is excited with any 
excitation shape y x . Nevertheless, when the disk rotates in water the structural mode +n is excited for 
an arbitrary number of exciters only if y x = n, which means that the excitation has to coincide with the 
structural mode shape in its shape and in its direction. When fi ^ n, the structural mode +n is only 
excited if the number of exciters accomplish Equation (14). 

3. Experimental Tests 

3.1. Test Rig Setup 

To study the dynamic behavior of a rotating disk structure (rotating in air and water) due to the 
excitation characteristic of the rotor-stator interaction experimentally, a rotating test rig setup has been 
developed. It consists of a disk connected to a variable speed motor. The disk rotates up to 8 Hz in air 
and in water. When the disk is rotating, it is excited with several excitation patterns created with four 
piezoelectric patches PI-876A12 attached in two different configurations (only one PZT is shown in 
Figure 5). The patches work in a range of -100 V to 250 V. The excitation signals are created with a 
NI-9263 Signal Generator (National instruments, Austin, TX, USA) using a Labview®-code. Four 
independent signals can be generated with this generator in a range of -10 V to 10 V. To work in the 
PZTs range an amplifier OEM 835 with gain 25 is used. This amplifier sends the signals to the PZTs 
and also to the Bruel & Kjaer Acquisition System (Naerum, Denmark). 

Figure 5. Rotating disk test rig. 



SLIP RING 




COMPUTER 




MOTOR 
CONTROLLER 











The response is measured with miniature and submergible Dytran 3006-A accelerometers (Dytran, 
Chatsworth, CA, USA), that are screwed on the disk. The added mass of the sensors and actuators is 
less than 1% of the mass of the disk, so it is checked that after installation of them, the vibration 
characteristics of the disk do not vary. 
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All the excitation and response signals of the rotating frame are transmitted to the stationary frame 
through a Michigan S10 slip ring system. This system is attached at the tip of the shaft. Ten 
independent circuit slip rings are used to transmit the signals. Since one channel of a sensor or exciter 
consist in a positive and a negative line, only four accelerometers and four patches can be used 
simultaneously. The negative terminals of the accelerometers are connected to one common point and 
the negative terminals of the patches are connected to another common point. The eight channels left 
are used to connect the four positve lines of the accelerometers and the four positive lines of the patches. 

The rotating speed of the motor (Mavilor MLV-072, Santa Perpetua de Mogoda, Spain) is 
controlled with a computer. The vibrations of the motor are isolated from the rest of the test rig 
through a silent block. A scheme of the test rig used in this study is shown in Figure 5. 

3.2. Rotating Disk and Piezoelectric Patches 

The disk is made of stainless steel with an external radius r out = 200 mm and a thickness of 
ho = 8 mm. The disk has a hole on its center in order to attach the shaft. The nomenclature used for the 
patches is P-X, where X is the angle related to the 0° direction in counterclockwise direction, when the 
disk is attached to the shaft and viewing the test rig from the top. The same nomenclature is used for 
the accelerometers, with A-X. Four accelerometers (A-0, A-90, A-135, A-180) are screwed on the disk 
to measure the response. Six patches PI-876A12 (61 mm x 35 mm) are glued on the disk (P-0, P-90, 
P-135, P-180, P-270, P-315) to make the excitation, but only four can be used simultaneously cause 
the limited number of channels in the slip ring system. In order to excite properly the studied mode 
shapes they will be used in two different configurations shown in Figure 6. Configuration 1 is (P-0, 
P-90, P-180, P-270) and configuration 2 is (P-0, P-135, P-180, P-315). The mass of the disk is 
approximately 7.6 kg. 



3.3. Proccedure 

The main objective of the tests is to study the response of the disk in resonance, when it is excited 
with different rotating excitation patterns that simulate the RSI excitation. This is done exciting the 
first several natural frequencies of the rotating disk in air and water with different excitation patterns y 
created with PZTs. The following procedure is applied: 



Figure 6. Rotating disk with installed PZTs and accelerometers. 




Sensors 2014, 14 



11930 



(1) Firstly, only one patch is used with a sweep signal that excites the first several natural 
frequencies of the disk. This is done for the disk rotating in air and the disk rotating 
in water. 

(2) It is checked that, when using the same excitation signal for two different patches the 
response of the contiguous accelerometer is different in terms of amplitude and phase, 
since the excitation depends on the mounting condition of the patch. Furthermore, for 
the same patch the relationship force/voltage changes within the excited frequency. 
Therefore patches have to be calibrated, in order to make a compensated excitation 
shape at one desired frequency. In this case, since the response of the disk is studied 
under resonance condition, the calibrated frequencies are the natural frequencies of the 
disk. For the first several natural frequencies, patches are calibrated to make a 
compensated excitation shape in angle and phase. 

(3) Once patches are calibrated, the disk is excited with several rotating excitation patterns 
at the studied natural frequencies. 

3.3.1. Excitation with One Patch 

Natural frequencies of the rotating disk in air and in water have to be determined. Since, the first 
several natural frequencies are the most relevant in the real case, cause they can be excited by RSI 
phenomena [14], this study is concerned in a frequency range of 0-1200 Hz. Therefore, for this disk a 
sweep signal from 0 to 1200 Hz is used to excite the first natural frequencies of the disk in air and in 
water. Such a signal can be described as: 

y = Asin(a)t(t)) for0<t<t end (15) 

When a patch works with this signal, it excites all the frequency band from 0 Hz to < —^ L Hz. 

< —^ L is selected as 1200 Hz in this case, co (sweep rate) has to be enough small (slow sweep) in order 

to have a good resolution in frequency when applying the FFT, without losing information. 

In Figure 7, the procedure to obtain the natural frequencies of the disk is shown for one resonance 
and one sensor. The time signal of the excitation P-0 (Figure 7a) shows a slow sweep excitation. The 
time signal of A-0 (Figure 7a) shows that a resonance occur at certain time. 

To obtain the frequency content of these signals, a Hanning Window of 4 s (resolution 0.25 Hz) is 
applied on the time signals. Since this window is shorter than the total length of the time signals, it is 
translated 0.2 s (5% of the window length) every average. In each average, the FFT is applied in both 
signals and superposed to the other averages with the maximum hold method, which considers only the 
maximum value for each frequency. In this way the frequency content of both signals is obtained 
(Figure 7b). Using both response (A-0) and excitation (P-0) signals, the frequency response is obtained 
(Figure 7c). Natural frequencies are detected in precision looking at the peaks of the frequency 
response function (FRF). The corresponding mode shapes are obtained analyzing the relative phase of 
the four accelerometers (Figure 6) and contrasting with the analytical model. In this way, natural 
frequencies and mode shapes are determined for the still and rotating disk in air and in water. 
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Figure 7. Determination of the natural frequencies of the rotating disk, (a) Time signals; 
(b) Signals after FFT; (c) FRF Amplitude and Phase. 
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3.3.2. Calibration of the Patches 



A closer view from the disk with the installed piezoelectric patches can be seen in Figure 6. The 
nomenclature used for the patches is P-X, where X is the angle related to the 0° direction in 
counterclockwise direction, when the disk is attached to the shaft and viewing the test rig from the top. 
The same nomenclature is used for the accelerometers, with A-X. 

The relation force/voltage characteristic for the piezoelectric patches changes for each patch (since 
it depends on the mounting condition of the patch) and with the signal frequency. To make that patches 
work with the same amplitude (in force) and with the desired phase to each other, they have been 
previously calibrated at the natural frequencies studied. Here is explained how the calibration is 
performed for the mode n = +2 (disk rotating in air) and configuration 1 (Figure 6). For other modes 
and configurations the procedure is equivalent. 
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First, only patch P-0 is used at one natural frequency with a peak value of 75 V. (XA-o(p-o))fn 
(amplitude of the vibration of accelerometer A-0 due to an excitation with P-0 at the natural frequency 
n) and ((XA-o(P-oyp-o)fh (angle between the signal of A-0 and the signal of P-0 due to an excitation with 
P-0 at the natural frequency n) are measured. 

When using another patch at the same natural frequency, (XA-i(p-i))fo and (aA-i(p-iyp-o) & are measured 
(i is 90°, 180° and 270° in this case). The amplitude of P-i is changed in order to accomplish 
(XA-i(p-i))f n =(XA-o(p-o))fh for each i. Also a phase shift between signal P-i and signal P-0 is introduced to 
make that (aA-i(p-i)/p-o)fn = ( a A-o(p-oyp-o)fn- In this case the signal of P-0 is acquired as a reference, but is 
not really used to excite the patch P-0. In the specified case, the calibration of the patches has been 
done adjusting the signals of the patches, to accomplish: 

(XA-0(P-0))fn = (Xa-90(P-90))& = (XA-180(P-180))fh = (XA-90(P-270))fn 

(16) 

(ctA-0(P-0)/P-o)fn - (ctA-90(P-90)/P-o)fn _ (ciA-180(P-180)/P-o)fn ~ (<XA-90(P-270)/P-o)fh 

The accomplishment of Equation (16) for each f n (f n are the first natural frequencies of the disk), 
guarantee that patches are properly calibrated in amplitude and phase (Figure 8b). After patches are 
calibrated, a phase shift between them (apart from the phase shift introduced for the calibration) can be 
introduced to make the desired excitation pattern. 

Figure 8. Polar plot of the sensors, (a) Before calibration of PZTs; (b) After calibration 
of PZTS. 




3.3.3. Rotating Excitation at the First Natural Frequencies 

After the first several natural frequencies are obtained and patches are calibrated to make a 
compensated excitation shape in these frequencies, the disk is excited with several rotating excitation 
patterns that simulate the RSI that occurs in real turbomachinery. This is performed for the disk 
rotating in air and for the disk rotating in water. The excitation patterns that are created with four 
patches installed in the two different configurations shown in Figure 6 are represented in Figure 9. The 
response is measured with the accelerometers screwed on the disk. 
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As Figure 9 shows, with four patches attached at 90° it is not possible to define the rotating 
direction of the excitation for the modes n = ±2 and n = +4, changing the phases between exciters. For 
n = +3 this direction is decided changing the phase of the patches (Figure 9). For n = +2 another 
configuration is tested (P-0, P-135, P-180 and P-315), that defines the rotating direction. With this 
configuration is also not possible to define the direction of the excitation for n = +4. 

4. Results and Discussion 

4.1. Natural Frequencies of the Disk 

The first several natural frequencies of the disk are determined for the case that the disk is rotating 
in air and in water. To determine these natural frequencies and mode shapes, two methods are used, the 
analytical method presented in Section 2 and the experimental method explained in Section 3. For the 
rotating case in air and water, the maximal velocity will be studied in this paper, i.e., when the disk is 
rotating at 8 Hz. 

4.1.1. Analytical 

To determine the natural frequencies in air for each structural mode +n, the density, thickness and 
averaged radius are necessary (Equation (5)). The parameter D* is determined using the methodology 
presented in [25]. In that study only the formulas for n < 3 are presented. Therefore, the selected mode 
to determine D* is n = +2. Once D* is obtained, the modes n = +3 and n = +4 are determined for the 
non rotating disk. These are the predicted natural frequencies below 1.2 kHz, which is the range that 
will be studied in this paper. According to the analytical model for low rotating speeds of the disk the 
rotation does not affect in Equation (5). Therefore the natural frequencies calculated for the non 
rotating case are also the natural frequencies for the disk rotating at maximal speed of 8 Hz. 

When the disk is submerged in water, apart from the parameters used in the case of air, the density 
of water, the distances to the rigid walls H up & Hd own and the velocity of the fluid with respect to the 
disk are needed (Equation (6)). To estimate the averaged values of Q. up and £2down, a CFD calculation 
has been performed. This calculation gives the velocity of the particles of fluid with respect the 
stationary frame. An averaged value for the rotating speed of the fluid particles with respect to the disk 
for the upper and lower field (Figure 1) is calculated. For the rotating speed of the disk of 8 Hz it is 
determined that £2 up = 4.97 Hz and f2 d own = 4.92 Hz. All the parameters used in Equation (6) (and also 
in Equation (5)) are presented in Table 1 . 



Table 1. Parameters used to determine the natural frequencies of the rotating disk in air 
and in water (£2 d i s k = 8 Hz) for each n. 



Longitudinal Parameters 


Densities 


Rotating Velocities 


Bending Stiffness 


(m) 


(Kg/m 3 ) 


(Hz) 


(N/m) 


r o H up Hdown h D 


Pd Pwater 




D* 


0.072 0.03 0.077 0.008 


7800 1000 


4.97 4.92 


253 
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4.1.2. Experimental 

From the experimental point of view, the disk is excited with only one patch for the stationary case 
and for the rotating case in air and in water. The way to obtain the natural frequencies of the disk 
below 1.2 kHz in all these cases is mentioned in Section 3. Figure 10a shows the natural frequency 
n = +2 of the disk in air for the stationary and for the rotating case and Figure 10b shows the same for 
the case that the disk is rotating in water. 

It is seen that for the case in air only a slight increase in the natural frequency is produced when the 
disk is rotating in air. This increase represents less than 1% for the studied modes. The phase shift 
between sensors (Figure 10a) shows that the corresponding mode shape is a standing wave (viewed 
from the rotating frame) with all the points moving in phase (0 rad) or in counterphase (71 rad) to each 
other. For the case that the disk rotates in water, an important difference is observed. In this case two 
peaks are detected (from the rotating frame), which correspond to the modes n = 2 and n = -2. These 
are two travelling waves travelling in the opposite direction as the analytical model predicts, and the 
phase of accelerometers confirms (Figure 10b). According to Equation (3), the phase shift between two 
points with a spatial phase shift of 0, have a phase difference of ±n0, depending on the mode shape 
and on the travelling direction. Therefore, for the shown accelerometers in Figure 10b (A-0 & A- 135), 
this phase shift has to be ±;t/2. The values of the natural frequencies of the disk below 1.2 kHz in all 
the cases for analytical and experimental methods are shown in Table 2. 

Figure 10. (a) Natural frequencies n = +2 of the rotating disk in air at Hdisk = 0 Hz and 
£^disk = 8 Hz; (b) Natural frequencies n = +2 of the disk rotating in water at Q.dkk = 0 Hz and 
^disk = 8 Hz. 
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Differences of less than 5% are obtained when comparing both methods. After the first several 
natural frequencies of the disk are obtained, patches are calibrated to make the desired excitation 
patterns in these natural frequencies. 



Table 2. Natural frequencies (Hz) of the stationary and rotating disk in air and water. 
Analytical and experimental values. 









Air 






Water 








Analytic 


Experim. 


Error 


Analytic 


Experim. 


Error 




n = +2 


260.03 


257.75 


0.9% 


151.37 


150.50 


0.6% 


Q disk = 0 Hz 


n = ±3 


585.07 


588.25 


0.5% 


393.20 


383.25 


2.5% 




n = +4 


1040.13 


1031.50 


0.8% 


750.85 


733.25 


2.3% 




n = +2 


260.03 


259.75 


0.1% 


145.33 
157.14 


143.25 
156.00 


1.4% 
0.7% 


f^disk = 8 Hz 


n = ±3 


585.07 


590.25 


0.9% 


386.23 
399.62 


374.00 
390.00 


3.1% 
2.4% 




n = +4 


1040.13 


1033.25 


0.7% 


743.25 
757.89 


708.75 
729.50 


4.6% 
3.7% 



4.2. Dynamic Behavior of the Rotating Disk in Air Due to an RSI 

The RSI is the superposition of excitations that change in frequency when increasing the rotating 
speed of the rotor (term ZZ 0 fi rot of Equation (8)) and with an excitation moving with respect the 
rotating disk (term kZ b — IZ 0 of Equation (8)). In this paper, the dynamic behavior of the disk when 
the excited frequency coincides with a natural frequency of the disk is studied. A slow sweep signal 
around the resonance with the patterns specified in Figure 9 is applied on the disk. Amplitudes of 
vibration of the accelerometer A-0 (attached on the rotating disk) under resonance condition are 
studied for all the cases represented in Figure 9. 

4.2.1. Analytical 

Applying Equation (11) particularizing the excitation shape and the structural mode shape the 
amplitude of the disk under different resonance situations is determined. 

4.2.2. Experimental 

Once natural frequencies are determined and patches are calibrated, the disk is excited with a sweep 
signal passing through the resonance with the excitation patterns presented in Figure 9. As Figure 1 1 
shows, for n = +2 when the disk is excited with the pattern y = ±2 the resonance is amplified and 
when it is excited with y = +4 is eliminated (as predicted in Equation (12)). Note that for the 
symmetric-position of patches the excitation direction for y = +2 and y = +4 cannot be defined. For the 
same structural mode, the disk is excited with the non-symmetric position of the PZT. In this case the 
direction of excitation can be defined for y = +2. As Figure 11 shows, both excitations amplify the 
amplitude of resonance, as the mode shape is a standing wave on the disk. 
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Table 3. Amplification of the resonances (A Y /Ai_ patch) of the rotating disk in air 
(£2disk = 8 Hz) due to the different excitation patterns. Analytical, experimental and error. 



Excitation Shape (y) 





y = 2 


y = -2 


7 = 3 


y = -3 


y = 4 y = -4 




2 


2 








n = +2 


1.981 
1.0% 

4 

3.942 
1.5% 


1.971 
1.5% 


0 

0.003 


0 

0.004 


0 

0.003 




0 




2 


2 


0 


n = ±3 


0.003 




2.005 
0.3% 


2.008 
0.4% 


0.005 




0 




0 


0 


4 


n = +4 


0.017 




0.003 


0.001 


4.042 



1.1% 



The rest of resonance amplitudes divided by the amplitude of resonance due to one patch excitation 
(red line in Figure 11) are presented in Table 3 compared with the analytical results. A difference of 
less than 1.5% between the analytical model and experimental results is obtained for all the 
experimented cases. From the experimental results can be extracted that patches are feasible to excite a 
thick disk in air and that RSI excitation can be studied and simulated (or attenuated) with PZTs. It is 
found that for each mode n, only one resonance (from the rotating system) is produced and amplified 
when the excitation shape y coincides with the structural mode shape n, with no matter of the direction 
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of rotation of excitation. For the other excitation shapes y, with the positions of the exciters used, the 
amplitude of resonance is almost 0. 

4.3. Dynamic Behavior in Water Due to an RSI 

4.3.1. Analytical 

Applying Equation (14) particularizing the excitation shape and the structural mode shape the 
amplitude of the disk under different resonance situations is determined. 

4.3.2. Experimental 

Again the same experimental procedure explained in Section 3 is applied when the disk is rotating 
in water. Figure 12 shows the excitation of the structural modes n = +2 and n = -2 with different 
excitation patterns. 



Figure 12. Resonances around n = 2 and n = -2 (experimental) for the disk rotating in 
water (£2disk = 8 Hz). Different excitation patterns. 




125 135 145 155 155 175 185 

Frequency (Hz) 



In this case, when the disk is excited with one patch, the two resonances (n = 2 and n = -2) are 
detected from the rotating frame. For the excitation y = ±2 without defining the rotating direction of 
the excitation, both resonances are amplified. When y is different than +2, the amplitude of the 
resonance is almost 0 as predicted in Equation (14) for the position of exciters used. Only one 
resonance appears if the rotating direction of the excitation is defined. In the case of the non-symmetric 
position of patches this direction is defined (Figure 9) and it can be checked that to excite the structural 
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mode n = 2 the excitation y = 2 is necessary. For the same structural mode the excitation y = — 2 
eliminates the resonance. The same conclusion is achieved for the structural mode n = -2, which needs 
an excitation shape y = -2 to be amplified. The rest of resonance amplitudes compared to the case of 
one patch excitation are presented in Table 4 compared with the results of the analytical model. 

Table 4. Amplification of the resonances (A y /Ai_ P atch) of the rotating disk in water 
(£2disk = 8 Hz) due to the different excitation patterns. Analytical, experimental and error. 

Excitation Shape (y) 



o 



u 
3 



in 





y = 2 




y = -2 


y = 3 


y = -3 


y = 4 y = -4 




4 




0 








n = +2 


3.852 
3.7% 


A 

4 

3.983 

U.4 /C 


0.120 


0 

0.007 


0 

0.008 


0 

0.003 




o 




4 










0.033 




3.991 








n = -2 




4 

A 1 A1 

4.101 

3% 


U.zvo 


A 
U 

0.003 


a 
U 

0.004 


A 
U 

0.015 






0 




4 


0 


0 


n = +3 




0.003 




4.023 
0.6% 


0.092 


0.007 






0 




0 


4 


0 


n = -3 




0.002 




0.025 


3.965 
0.9% 


0.005 






0 




0 


0 


4 


n = +4 




0.005 




0.009 


0.007 


4.152 

3.8% 






0 




0 


0 


4 


n = -4 




0.003 




0.008 


0.017 


4.102 



2.6% 



A difference of less than 4% between methods is obtained for all the tested cases. From 
experimental results can be extracted that patches are feasible to excite a thick disk in water and that 
the RSI can be studied and simulated (or attenuated) with PZT. To notice is that the disk has two 
natural frequencies for each mode +n (viewed from the rotating system) when it rotates in water. These 
are detected when excited with one patch. The resonance is amplified only in case that the excitation 
shape y coincides with the structural mode n in magnitude and rotating direction. Both resonances are 
amplified when the rotating direction of the excitation shape is not defined. For other excitation 
patterns (n 4- j), with the used exciters, both resonances are eliminated. 
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5. Conclusions 

The dynamic behavior of a thick disk rotating in air and inside a casing filled with water has been 
analyzed experimentally and analytically. The disk has been excited with several rotating excitation 
patterns simulating the rotor-stator interaction (RSI) excitation. 

For an accurate analysis of the disk behavior both exciters and sensors measuring the response were 
located on the disk (rotating frame). For the excitation several PZT actuators attached to the disk were 
used. PZTs do not affect the mass of the disk and do not perturb the flow of water produced by the 
rotation of the disk inside the casing. PZT actuators have been used several times to excite thin rotating 
disks in air and from the stationary frame, but never to excite a thick disk submerged in water and 
confined inside a casing. Experiments presented have demonstrated that it is feasible to use PZT 
actuators in thick structures in air and submerged in water. This is interesting because PZTs could be 
used to determine the dynamic response of disk-like structures, such as turbomachinery impellers in 
actual operating conditions. 

Natural frequencies of the rotating disk in air and in water have been determined. The effect of 
rotation is very important when the disk rotates in water. Viewed from the rotating frame, even for low 
rotation speeds (8 Hz), for each natural frequency found in the non rotating case (one for each 
diametrical mode), two natural frequencies that correspond to two travelling waves rotating in the 
opposite direction appear. For the same rotating speed in air, only one natural frequency that 
correspond to a standing wave (viewed from the rotating frame) is observed. This behavior of the 
natural frequencies and mode shapes (when the disk rotates in air and in water) is determined 
experimentally and predicted by the analytical model used. 

The dynamic response of the rotating disk in water at one natural frequency depends on the 
excitation shape and also on the rotating direction of the excitation. When the excitation is fixed on the 
rotating frame, the response at the two natural frequencies of the corresponding diametrical mode is 
amplified. However when the excitation spins in the same direction of the disk only the lower natural 
frequency is excited, which corresponds to the travelling wave travelling in the same direction as the 
disk. When the excitation spins in the opposite direction only the higher natural frequency is excited, 
which corresponds to the travelling wave travelling in the opposite direction. If the excitation shape 
does not coincide with the diametrical mode considered, the response at both natural frequencies is 
almost zero (for the studied configuration of patches). When the disk rotates in air, only the excitation 
shape affects the dynamic response. If the excitation shape coincides with a diametrical mode, the 
response of the disk is amplified at the corresponding natural frequency, with no matter of the rotating 
direction of the excitation (rotating with the disk, counterwise or standing). These results are obtained 
experimentally and using the analytical model. 
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